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1 Introduction. 



Unconventional properties of the oxides like cuprates (Y Ba2Cu^O^ +X) La2- x Sr x CuO& y La2CuO^ + d,.. 
.), manganites (Lai- x Sr x MnC>3,...), nickellates (Lo^iVzQi+d,...), bismuthates ((K, Ba)BiO^) 
including systems with the high-T c superconductivity and colossal magnetoresistance reflect 
a result of a response of the system to the nonisovalent substitution or charge disorder 
that stabilizes the intermediate valence phases providing the most effective screening of the 
charge inhomogeneity. These phases in oxides may involve novel molecular cluster configu- 
rations like the Jahn- Teller sp-center IJ with anomalously high local polarizability. Their 
appearance may be a result of the disproportionation reaction like 

M + M = M + +M~, 

where M is a basic metal-oxygen center (CuOf', NiOf~, Mn0 9 f , BiOl~ , respectively), the 
M centers are the corresponding polar (hole and electron) centers. All these oxides are the 
so called charge transfer semiconductors, where a fundamental absorption band is determined 
by the charge transfer s — > p transition from the nondegenerate and predominantly metallic 
even s-state to the low lying degenerate and predominantly oxygen odd p-state. An ionization 
or a hole doping for the sp-like M-center may be accompanied by the quasi- degeneracy effect: 
a hole may be localized either at the predominantly metallic (s) or at the predominantly 
oxygen (p) molecular orbital with competition of two configurations, s 2 and sp, respectively. 
As a result, we have for the hole center M~ a ground manifold of the terms with different 
parity, spin multiplicity and orbital degeneracy providing a multi-mode behavior of such 
centers. First of all, these are unstable with respect to the pseudo- Jahn- Teller effect [[IJ [2|. |3| 
with active even and odd local displacement modes of different symmetry and formation of 
the appropriate vibronic states |fj]. The hole pseudo- Jahn- Teller (PJT) center with its high 
polarizability can be a center of an effective local pairing with a formation of the local singlet 
(e.g. cuprates) or triplet (e.g. manganites) boson as two electrons paired in molecular shell. 
Thus, we come to the electron PJT center: M~ = M + + (local boson). A local boson can 
correspond to the completely filled molecular shell with the 1 A lg symmetry (e.g. b\ g for the 
cuprates and nickellites and 6s 2 for the bismuthates) or to the half-filled molecular shell e 2 
with the 3 A 2g symmetry in manganites. Two polar centers may have a similar structure of 
the ground manifold. It is worth to note that the PJT nature of the polar centers results in 
a strong vibronic reduction for a probability of their recombination. As a whole the oxides 
under such consideration may be called as the strongly correlated PJT oxides. 

In particular, a model approach developed in papers JI], ||, [| considers the C11O4 cluster 
based copper oxides as systems unstable with respect to the disproportionation reaction 

2CuO\- -> [CuOl-]pj T +[CuOl~]p JT 

with the creation of the system of the polar hole (h) CuO\~ or electron (e) CuO\~ pseudo- 
Jahn- Teller centers. These centers are distinguished by the local boson or by two electrons 
paired in the completely filled molecular orbital of the CiiC^-cluster. The new phase can 
be considered as a system of the local bosons moving in the lattice of the /i-centers or as a 
generalized quantum lattice bose-gas. An origin and anomalous properties of the /i-centers 
are connected with a near degeneracy of the molecular terms l A\ g (Zhang- Rice singlet), 
1 E U , 3 E U for the configurations b\ g and bi g e u , respectively, that can create conditions for 
the pseudo- Jahn- Teller effect with active local displacements modes of the Q eu , Qb lg and 
Qb 2g types. In general, the PJT effect in the 1 Ai g , 1,3 E u manifold leads to the formation of 
the four-well adiabatic potential of two symmetry types: B\ g or B 2g with the nonzero local 
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displacements of the Q eu , Qb lg or Q eu , Qft 2g -types, respectively (Fig.l). Unusual properties 
of the 1 A ig , 1,3 E U manifold with the terms distinguished by the spin multiplicity, parity 
and orbital degeneracy provide unconventional behaviour for the ^-centers with an active 
interplay of various modes. The hole PJT center with its high polarizability can be a center 
of an effective local pairing. 

A most direct and convincing observation of the hole singlet-triplet PJT centers is made 
recently H by the NQR method in La^Cu^Li^O^. The authors revealed the spin singlet 
ground state [S = 0) and the low lying spin triplet state (S = 1) with the singlet-triplet sep- 
aration Ast = 0.13ey comparable with the well known estimates for the nearest neighbours 
exchange interactions for the parent cuprates. They found the anomalously weak temper- 
ature dependence of the relaxation rate at low temperatures that evidences an occurrence 
of the spinless multiplet structure in the Cw0 4 cluster ground state. They revealed an ap- 
preciable spin contribution to the low temperature relaxation indicating the simultaneous 
occurrence of the ground state multiplet structure, the sufficiently low singlet-triplet sepa- 
ration and the intrinsic singlet-triplet mixing. These features and the observed relaxation 
inequivalence of the various Cu sites are quite natural for the PJT centers in the conditions 
of the static JT-effect. 

The PJT hole centers like CuO\~ with the singlet-triplet quasi-degeneracy within a 
ground state have been observed by ESR-spectroscopy in LaSr Al\_ x Cu x O \ which is isostruc- 
tural to La2- x Sr x Cu04 |]]. 

A description of the multimode PJT centers system is very complicated. The simplest 
model Hamiltonians correspond to the quantum lattice bose-gas Hamiltonian for the charge 
subsystem M, M and the Hamiltonian of the cooperative JT-problem for the vibronic sub- 
system Q. Below we'll present a simplest model approach to the spin subsystem of the PJT 
centers phase. 

Unconventional properties of the spin subsystem in the PJT centers phase as for the 
generalized singlet-triplet magnetic system with well developed charge, spin and structure 
fluctuations are connected with the induced spin fluctuations. The main source of them is 
determined by a dependence of the singlet-triplet separation A ST on the local boson density. 
The local boson movement is accompanied by the fluctuations of the Ast value and appro- 
priate spin fluctuations. In general, for a spectral region of the spin system transparency (or 
far from the main spin excitations) the temporal dependence of the induced spin fluctuations 
may be entirely determined by the charge fluctuations. This enables to explain many un- 
usual spectral and momentum (u>, q), temperature and concentration (x, T) dependences for 
different effects determined by the spin correlation functions including an inelastic neutron 
scattering (INS) and the spin-lattice relaxation. Moreover, it indicates a principal possibility 
to make use of the magnetic methods (INS, NMR, spin susceptibility measurements) for the 
revealing and studying the non-magnetic (charge and vibronic) excitations. 

The paper is organized as follows. In Section II we consider the mean field approximation 
(MFA) for the simplified singlet-triplet Hamiltonian of the PJT centers system with total 
neglect the intermode coupling. In Section III we suggest a model approach to account of 
the spin-charge correlations within a so called induced spin fluctuations model. Some static 
and dynamical spin properties for various copper oxides are discussed in Section IV in a 
framework of the singlet-triplet model. 
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2 The singlet-triplet model of the polar PJT centers 
phase. Mean field analysis. 

To maximally simplify the spin subsystem description within the PJT centers phase we'll 
restrict ourselves only to the singlet and triplet 1 E U , 3 E U terms for the &i s e u -configuration (see 
Fig. 2) neglecting the orbital degeneracy and taking into account the isotropic spin exchange 
interaction between b\ g and e u holes of the neighbouring C11O4 centers. Thus, we come to 
the Hamiltonian of the two-component spin liquid 

1 1 1 

H ST = g 2 h^eS^Sli^i + ^ H J&i 9 &i 9 (u)siiSy + 2 H h lg e u {ij)siiS 2j + 
i ij ij 



+ 2 S ^^(uO^sij + 2 H J eiieil (u)s2iS2i, (1) 

where the s\j, s*2i are the spins localized on the predominantly copper 6 lg or the pure oxygen 
e u orbitals, respectively. First term in (1) describes the intra-center c/p-exchange and the rest 
ones describe the inter-center dd, dp, pjo-exchange, respectively. Introducing the summary 
spin of z-th cluster 

Si = Sn + S2i 

and the spin polarization operator 

Vi = su - s 2i 

with the non-zero matrix element like (00 1 V z | 10) = 1 responsible for the singlet-triplet mixing 
one can rewrite (1) as 



Here we used the new exchange integrals 

l\f = \{KhM + Wtf) + 2h lg eM)}, (3) 

f = 4 s) = ^w«)-uw)}. ( 4 ) 

4f = tUW^') + Wu) - 2h lg e u (ij)} (5) 



and also introduced the singlet-triplet separation Ast = h lg e u (H)- The first term in (2) 
explicitly takes account of the intra-center dp-exchange resulting in the singlet-triplet sep- 
aration for each cluster, the second one describes the conventional Heisenberg exchange of 
triplets leading to the formation of spin waves (SW). The rest terms are non-trivial and 
correspond to the important effects of spin polarization exchange with possible formation of 
magnetic excitons (spin multiplicity waves) and their hybridization with SW. 

Numerical values and signs of the exchange parameters jC 1 ' 2 - 3 ' 4 ) are determined as a result 
of a competition among the b\ g — bi g , b\ g — e u and e u — e u exchange interactions. Take notice 
of different signs of the b\ g — e u contribution to and 1^ and also of the b lg — b± g and 
e u — e u contributions to 1^ = 1^ . 
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Contrary to the well known weakly covalent and weakly correlated magnetic oxides like 
ferrites, the exchange interaction for two nearest neighboring PJT centers in the strongly 
covalent cuprates || includes the unusually large ferromagnetic Hund (for different orbitals) 
or antiferromagnetic correlational (for similar orbitals) contributions of the intra-atomic Oi v - 
exchange within the common oxygen ion. Besides, we should deal with a relatively large 
ferromagnetic contribution of the Heisenberg Cu3d — 02p exchange. As a whole, a detailed 
microscopic analysis of the exchange interactions for the PJT centers requires a special 
consideration, especially, in conditions of their comparable large values / ~ Ast ~ 0.1 eV. 

Notice that some kinds of the singlet-triplet Hamiltonians have been discussed to describe 
the magnetic behaviour of some rare-earth magnetic systems with Pr 3+ , Tb 3+ , Tm 3+ ions [fT0| , 
[TTJ. Recently, the singlet-triplet magnet model has been suggested for La 2 A r «0 4 compound 

Above we neglect the nonmagnetic Zhang- Rice singlet 1 Ai g which inclusion in the H$t 
leads to an essential complication with a renormalization of the exchange parameters. Note 
only an appearance of the supplementary contribution like 



J2 

ij 



as a result of the inter-center electrostatic interactions U with account of the different electron 
density distribution for singlet 1 A\ g and triplet 3 E U states. This term leads to the possibility 
of the singlet-triplet ordering (like the charge one) without any spin-exchange interactions. 



The Hamiltonian (2) has the SO (A) group symmetry JT3[ 

[Si, Sj] = ieijkSk, [Si, Vj] = iEijkVk, [Vi,Vj] = zeijkSk, (6) 

an algebra of which has two invariant operators (Casimir's operators) 

Cf ° (4) = S 2 + V 2 = 31, C s 2 ° {4) = S-V = 0. (7) 

These relations provide a kinematic constraint for S and V modes and automatically follow 
from the ordinary spin algebra. 

Below, we'll treat the Hamiltonian (2) in the framework of the modified MFA approach 
within the nearest neighbours approximation with two sublattices labelled by additional 
7 = A, B indices. Then Hamiltonian (2) can be rewritten as 

HmFA = 5^ ^i-y ~ 2 5-^ ^sia-y(Si a )j ~ ~ h v icr( (^ia) 7+ (8) 
27 iay iay 

-f- ^ ' hgiaySicfy ~\- ^ h v i a ~jVi a ~j (9) 



ict-y icey 



with molecular fields 



hsia, = + l8 } (V ja )j), (10) 

hvurr = E(4 2) (^)y + 4 4) <^>y)- ( n ) 
3,1 

Note, that the MFA treatment of the Hamiltonians (1) and (2) is different in what 
concerns the taking account of the first term describing the intra-center exchange interaction. 
Contrary to the usual MFA, a summary spin S'M^-representation used in the Hamiltonian (2) 
provides an explicit account of the intra-center exchange. Within the usual MFA approach 
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— * — * — * — * 

for the Hamiltonian (1) the averages (S) and (V) (often denoted as M and L) correspond 
to the intra-center "ferromagnetic" and "antiferromagnetic" vectors for the PJT center, 
respectively, with kinematic constraint 



(Sf + (V) 2 = As\ = Ash (S)-(V) = 

instead of the operator constraint (7) which in the appropriate averages is equivalent to 

(S 2 ) + (V 2 ) =4si(si + l) =3, 0-V) = O. 

So, the singlet-triplet model within the modified MFA approach permits more correct account 
of the quantum effects essential at low spin magnitudes s± t 2 = \- 

An occurrence of two vector order parameters even with kinematic constraint (7) makes 
an analysis extremely complicated. The relatively simple collinear spin configurations like 
those that depicted in Fig.3(b,c), could be obtained only for four particular cases: 1) 

7(1,2,3,4) < 0> = = = 0> 2) J( l,4) > 0> /( 2,3) < 0> Q AA = ^ = = ^ 

3) I&W* > 0, Q aa = Qbb = 0, Qab = tt, 4) < 0, 1^ > 0, Qaa = Qbb = n, 

Qab = 0, as direct analogues of the ferro- and antiferromagnetic ordering. Take notice of 
the corresponding SS- and VF-configurations could be obtained as the usual MFA solutions 
of the Hamiltonian (1) only at T = when the MFA constraint (S) ■ (V) = is fulfilled both 
for the SS ((V) = 0) and the VV ((S) = 0) phase. 

In general, we have to deal with a frustration of the SS, VV, SV type exchange inter- 
actions that implies an appearance of the two-sublattice multiparametric ((S), (V),Qaa = 
Qbb, Qab) angular S, V spin configurations schematically depicted in Fig. 3(a). Their occur- 
rence extends substantially a set of admissible spin states and results in effective involving of 
the spin system to the strong inter-mode coupling and hybridization that raises its role in the 
general optimization procedure, especially, with taking account of the large inhomogeneity. 

Leaving the detailed MFA analysis of the phase diagram for the singlet-triplet Hamil- 
tonian (2) extremely complicated for the angular modes to the separate consideration we'll 
consider at first a specific case of the collinear S, V structures. 

In the nearest neighbours approximation for the collinear spin structures the free energy 
per site is 

F=Tj = -- a E (h szl (S z ) 1 + h VZJ {V z ) l ) + ^- 

iv 4 y=A,B z 

E ln{2(cosh(^a 7 )+exp(-/3^)cosh(/3^ Z7 ))}. (12) 

Here a 7 = {A 2 + 4/i 2 Z7 } 5 , (5 — ^ denotes the inverse of the temperature. 

The set of self-consistent equations for the order parameters minimizing F is 



exp (— sinh {(3h s 
cosh (f ai T ) + exp (— cosh (Ph szl , 



iq\ CA ^l r J 2' 011111 \h" b sz~t) / 1Q s 
\O z ) - — -S — " — — -, [i-6) 



For T = this system gets a form 



2^sinh(fa 7 ) 

(y z ) = ^ V2 " . (14) 

7 cosh (| <Xy) + exp (— cosh {(3h szl ) 



_9^£7 I-j > \U I _ A 

^ a 7 ' 2 "7 ^ \ u sz-y\ 2 
j 2~Q!-y *\ | hsz'y | ~2~ 
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^ ^ ~ I -a > l/i I - - ' 

Obviously, these equations have the following non-trivial solutions: 

1) SS-phase: ({S Z ) A = Vl ^ 0, (^) B = % ^ 0, (\/) A = (V Z ) B = 0), where ^ = -r/ 2 = 
±1 > 0) or 771 = 7/2 = ±1 < 0). This solution is realized when 

A<2|/i S27 |-a 7 (7 = A, B) (15) 

or with account of the appropriate order parameter values 

A + {A 2 + (22/ (2) ) 2 }" <2z\I^\. (16) 

(Here and below z is the nearest neighbours number, i~( 1,2 ' 4 ) is an effective nearest neighbours 
exchange integral). This relation means that one of the pure triplet states |1 — 1) or |11) has 
the minimal energy at both sublattices A, B. 

2) VV-phase: ({V Z ) A = m ^ 0, (V Z ) B = r] 2 ^ 0, (S Z ) A = (S Z ) B = 0). By similar way we 
obtain 771 = T) 2 = v U (4) < 0); Vi = -V2 = V U (4) > 0); 

This kind of solutions is realized if 

A>2|/i S27 |-a 7 (7 = A, B), (18) 
that corresponds to nonequality 

\I^\{2z\I^\ + A} > \I^\{(2zI^) 2 - A 2 }\ (19) 

It means that on each site A or B the state minimizing the energy is a quantum mixture of 
the 1 00) and 1 10) states. It is worth to note that the situation widely explored in theoretical 
investigations (well isolated singlet) is the particular version of the ^F-phase. Only when the 
singlet-triplet separation is large enough compared with an effective magnetic polarisation 
exchange (A > 2zl^) the spin subsystem of the new phase of the polar PJT centers has 
the singlet ground state. 

— * — * 

Note again, that within the MFA approach the S and V vectors correspond to the "fer- 
romagnetic" and " antiferromagnetic" vectors for the C-u0 4 -center, respectively. Fig. 4 shows 
qualitatively an orientation of the copper and oxygen spins for different S, V arrangements. 

Above we have considered the various SS and VV modes at zero temperature. At T 7^ 
we have, in general, the nonzero values both for (S z ) and for (V z ) or, in other words, the 
phase mixing (SS — VV) occurres. The unconventional non-Brillouin temperature behaviour 
of the appropriate order parameters and spin susceptibilities for the obtained phases are 
presented in Fig. 5, 6 together with the simplified spin level schemes for the antiferromagnetic 
SS and VV ordering (exchange parameters i^ 1 ' 4 ), A$t = h lg e u (i'i) are chosen to be positive, 
exchange parameters f( 2 ' 3 ) are chosen to be negative). 

The Curie (Neel) temperature T c is defined as a temperature when all spin order param- 
eters tend to zero simultaneously. Expanding the system (13-14) on the corresponding small 
order parameters we get the pair of nonlinear equations 

2 exp (/?— ) sinh ((3— ){(2zl^) 2 - 4z 2 /W J< 4 >} = (3 + exp (J3A))x 
2 2 



{^(3 + exp (/3A)) ± 2zlW A ± 4^exp sinh (20) 

from solution of which we get Ti j2 , then T c = max{Ti, T 2 }. 

An existence of two spin order parameters leads to many unconventional properties of the 
cuprates. First of all, outline various possible types of magnetic arrangement in the system 
of the PJT polar centers. Along with the well studied phase of (anti)ferromagnetic ordering 
of the Neel type ((S z ) ^ 0, (V z ) = 0) but with the non-Brillouin temperature dependence 
the new phases can occur. Only when (S z ) = (V z ) = 0, (S 2 ) = 0, (V 2 ) = 3 (T = 0) all 
polar centers have the singlet ground state. Non-trivial phases are the solutions with the 
non-zero order parameter (V z ) when (S z ) = 0. These correspond to different spin density 
arrangements. 

Unusual low temperature dependences of the order parameters and susceptibilities may 
be revealed by the singlet-triplet system in the SS — VV crossover regime, when 

A pa 2\h szl \ — a 7 (j = A,B). 

An occurrence of two spin order parameters results in two types of the spin susceptibili- 
ties: the usual one 

d(S) 

and a new spin-polarisation susceptibility 

Xv ~ — — 
dh 

with a principally different temperature behavior presented in Figs 5,6. Note that the MFA 
predicts new spin-polarization susceptibility xv should turn into zero within paramagnetic 
region as a consequence of the zero's values of the corresponding molecular fields and an 
absence of the singlet-triplet mixing. So, experimental reveal of the nonzero xv w iU indicate 
a relative role of the W-correlations. 

Let note that the taking account of the singlet-triplet mixing or the V-terms results in 
a suppression of the "S'-type" magnetism with a peculiar quantum "nearly paramagnetic" 
behavior of the system. This reflects both in a behavior of the iS-order parameter and in a 
tendency to coincidence of the transversal and longitudinal susceptibilities. 

The partial copper and oxygen spin polarizations within the PJT center and the ap- 
propriate susceptibilities can have unconventional behavior for the considered phases de- 
termined by the antiferromagnetic SS and VV exchange interactions. So, with increasing 
the temperature for the S'S'-phase we have a crossover from the ferromagnetic intra-center 
ordering ((V z ) ~ 0) to the state with intensive quantum suppression of the copper polar- 
ization ((V z ) ~ —(S z )). In the UV-phase a similar suppression occurs beginning from the 
antiferromagnetic intra-center ordering. Both phases reveal unconventional behavior of the 
copper Xd = \{.Xs + Xv) and oxygen Xv = \(.Xs — Xv) susceptibilities. The close values 
and the same signs of the xs and xv result in a large magnitude of the former and a near 
cancellation of the latter. Naturally, in paramagnetic state both susceptibilities coincide. 



3 Noncollinear phase. 

In considered model the principal feature of noncollinear spin configurations is their two- 
sublattice nature with combined ferro-aniferromagnet arrangement of vector order parameter 
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components (S) and (V") (see Fig. 3(a)) and wave vector Qaf = (tt, tt)- This fact distinguishes 
them from the usual incommensurate noncollinear spin structures. 

The analytical consideration of angular phase is cumbersome and hardly advisable there- 
fore we discuss only the numerical results. The temperature behaviour of longitudinal and 
transversal order parameters for noncollinear phase is presented in Fig. 7(a) (for the sake 
of convinience we have choosed z axis along (S Z ) A vector). In contrast to collinear phases, 
in considered situation both order parameters (S) and (V) are nonzero simultaneously at 
T = though the orthogonality condition 

(S-V) = 

is satisfied for each site. 

The further temperature evolution of spin arrangement resembles the well known orien- 
tational transition in antiferromagnets. In Fig. 7(b) we present the temperature dependence 
of Qaa = Qbb angles between (S) and (V) vectors on A, B sites, Qab between (S) A and 
(S) B vectors. From this plot one can see that with increasing the temperature the sub- 
lattice turning occurs when the spin parameters attempt to be collinear (S) A \\ (S) B and 
(Y) A II ( (explicit type of arrangement is determined by the signs of exchange interac- 
tions J* 1 ' 2 ' 4 )). Be ginning from the temperature T* on-site spin order parameters (S) and 
(V) line up along z axis that is accompanied by sharp decreasing to zero of appropriate 
transversal components. It should be noticed that for noncollinear phase a clear difference 
between longitudinal xl (in Fig.8 applied magnetic field H \\ z) and transversal xt ( m Fig. 8 
H _L z) susceptibilities disappears that determines nearly paramagnetic behaviour of system 
up to low temperatures and looks like spin-glass behaviour. 

As for collinear phases xh Xv turn zero a ^ Neel temperature. We pay attention to 
anomalous temperature behaviour of the spin polarization susceptibilities x\'~ L - At natural 
positive xJI'- 1 the spin polarization susceptibility x\ changes sign at T ~ Q.lSeV and Xv * s 
negative everywhere. 

In connexion of full energy optimization, an appearance of multiparametric angular spin 
configurations instead of the rigid collinear (anti)ferromagnetic structures essentially extends 
the possibilities of spin subsystem in cuprates , especially, in conditions of strong inhomo- 
geneity of intra- and intercenter exchange parameters and complicates an observation of 
magnetic arrangement in experiments. 



4 Induced spin fluctuations. 

Of course, the above considered singlet-triplet Hamiltonian within the mean field approxi- 
mation gives an oversimplified picture of the real situation in cuprates and should be applied 
to the real systems with some caution. The developed singlet-triplet model corresponds 
actually to the weak inter-mode coupling regime for the cuprates with relatively independent 
charge, spin and vibronic subsystems. Within this model approach an inter-mode coupling 
can be reduced to the renormalization of the parameters such as the transfer integral, the 
exchange integrals etc. So, in full analogy with vibronic reduction [14|] we could introduce a 



spin reduction of the local spinless boson transfer integral resulting in more or less effective 
suppression of the critical temperature T c . Note that within a weak inter-mode coupling 
regime this suppression is determined only by the dependence of the spin state of the PJT 
center on the boson density and is slightly dependent on the character of spin-spin corre- 
lations. In general, a magnitude of the spin reduction of the local boson transfer integral 
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and its dependence on different factors such as concentration can essentially influence the 
superconducting properties of the cuprates. 

In general, for the cuprates we should account for the strong intermode spin-charge- 
vibronic coupling and an occurence of the strong charge and structure fluctuations resulting 
in a suppression of the long range spin order. Actually, we deal with complicated multi- 
component strongly correlated quantum bose-liquid with the hybrid charge-spin-local struc- 
ture fluctuations providing the most effective screening of the charge inhomogeneity. Nev- 
ertheless, the cuprates can exhibit the various short range spin fluctuations revealed in the 
magnetic inelastic neutron scattering and the spin relaxation phenomena. 

Among the various intermode coupling effects we mention shortly an influence of the 
vibronic PJT-modes on the exchange parameters and in more details a dependence of the 
singlet-triplet separation A$t on the boson density. 

Fig. 2 clearly indicates the strong dependence of the exchange parameters on the vibronic 
modes that can be expressed within a pseudospin formalism, for example, as 

I(mn) = I + Ii(a z (m) + a z (n)) + I 2 a z (m)a z (n) + I orb i z (m)i z (n) 

where l(m) is the Izing like orbital moment occuring for the e u -hole. Parameters and I orb 
are determined by the e u -hole contribution. Thus, the spin fluctuations should be strongly 
coupled not only to the vibronic or local structure quadrupole fluctuations but to pair local 
orbital current fluctuations. 

Quantitative account of the strong intermode coupling in cuprates is an extremely com- 
plicated challenging problem which needs in some model simplifications. Below, within the 
so called induced spin fluctuations approximation we'll consider in details the role of the 
well developed static and dynamic charge fluctuations. Within the local spinless bosons 
model and weak inter-mode coupling regime their effect on the spin subsystem is mainly 
determined by the dependence of the magnitude of the singlet-triplet separation A ST on the 
boson density distribution 

A ST (m) = A ST (m) + £ A(mn)SN(n), (21) 

n 

where A 5r (m) is the magnitude of the singlet-triplet separation for the m site in the absence 
of the charge fluctuations, 5N(n) = N(n) — (N(n)) is a boson density (number) fluctuation. 
The second term in (21) with A{mn) is a contribution to the singlet-triplet separation for 
the m site determined by the charge fluctuations on the n site. 

Parameters A{mn) at m = n describe the local intra-center correlation effects and those 
atm^n describe the nonlocal crystalline field effects. Within the above mentioned simpli- 
fied two-terms ( 1,3 E U ) model the singlet-triplet separation 

A st = E( 3 E U ) - E{ l E u ) 

has to be slightly dependent on the crystalline field effect and an appropriate non-local con- 
tribution to A,st is rather small. Within a real three-terms ( 1 Ai s , 1 ' 3 E u ) model a ground 
singlet is a vibronic mixture of the singlet terms 1 A ig and 1 E U , so the singlet-triplet sepa- 
ration A S t should be strongly dependent on the crystalline field effects. The fluctuational 
contribution in the crystalline field parameters can be represented in terms of the partial 
crystalline field parameters 

AB kq = J2h q (R mn )SN(n), (22) 

n 
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that, in their turn, can be reliably evaluated, for instance, within a point charge model 



2 e 2 

bkq(Rmn) = pfc+l (^mn) ■ (23) 
mn 

The fluctuations of the singlet-triplet separation result in the appropriate spin fluctua- 
tions for which in the framework of the mean field approximation we'll obtain a relatively 
simple linear coupling 

6S m (t) = E(Ptnnh m (t)+ P l n h v m (t))6N n (t), (24) 

n 

where h s,v is a molecular field of the S or V type, respectively. The known expressions for 
the molecular fields (10)-(11) enables us to couple the spin fluctuations 5S m (t) with the S 
and V operators 

6S m (t) = J2( p m nS m (t) + P^ n V m (t))6N n (t), (25) 

n 

where the quantities P^ represent actually the spin-charge susceptibilities. Molecular fields 
and spins in the right hand side of the expressions (24) and (25) have been taken at the 
average boson density without an account of charge fluctuations. 

The spin-charge susceptibilities pf£ n and P m 'n, in general, can vary within sufficiently wide 
range. In this connection, let outline the specific case of the so called optimized systems with 
identical structure of the low lying terms for the electron and hole PJT centers [|l4j]. Within 
this "rigid ground multiplet" approximation A(mn) = 0, P^ n = and SS m (t) = 0. 

The spin fluctuations SS m (t) determined by the charge fluctuations in accordance with 
(25) will be further called as the induced spin fluctuations. These fluctuations reflect 
the charge-spin mode hybridization that determines their unconventional ^-dependence even 
for the simple charge-spin modes with the spin structure like spin density wave (SDW) 

S m {t) = S (t) cos (q sp R m ), (26) 

specified by the wave vector ±q sp , and, in general, for the incommensurate charge structure 
like charge density wave (CDW) 

SN m (t) = 5N (t) cos (q ch R m + a), (27) 

where q c h is a specific wave vector and a is a relative CDW-SDW phase shift. 
In common case, we can convert (25) to the expression 

= jf E Pi^Q(t)8N^(t) (28) 
Q 

with certain relation for the appropriate q components. (Here and further we restrict ours- 
ervels with the " S'-type" contribution. The corresponding account of the " V-type" contri- 
bution is trivial.) 

Below, we assume this relation is valid not only for classical quantities S, SN, SS but 
for the corresponding operators. The appropriate correlation functions for the induced spin 
fluctuations can be represented as 

(6fa0)6§-&)) = ± £ P^ 0I x 0^(O)§^(t)SN^O)SN_^{t)) t (29) 

Q,Q' 
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where the hybrid spin-charge correlation functions are presented. An approximation of 
weakly interacting charge and spin subsystems for (29) implies a breaking apart the hybrid 
correlation functions, i.e. 

(4-q(°)^- ? - + Q'W^q(°)^-Q'(*)) « (4-q( )^Q'W)(^q(°)^-Q'W)- (3°) 

An appearance of the hybrid correlation functions in (29) directly points to the possibility 
to detect the non-magnetic charge dynamic fluctuations by means of the magnetic methods 
(inelastic neutron scattering, spin resonance, nuclear resonance, etc.). Let note an appear- 
ance of the critical behavior of the charge fluctuations near the charge ordering temperature 
Tco within the local bosons system (see Fig.9). This unconventional behavior of the charge 
system will be revealed in unusual temperature behavior of the induced spin fluctuations 
and in the appropriate spin properties. In particular, note the pseudo-gap behavior of the 
induced spin fluctuations below the temperature T co of the charge ordering in the boson 
system. 

Above we restricted ourselves with a case of the single C-uC^-plane cuprates like La2- x Sr x CuO± 
with one CwCVcenter per unit cell. For the "bi-plane" cuprates like Y BaiCu?P§+ x it is nec- 
cessary to take account of two CwCVcenters per unit cell. Then modified expression (29) 
for the induced spin fluctuations takes the form 

(^0)^t)) = i^ £ e^e-^'^'P s (Q)P s (-Q')x 

Q,Q' u,v'=l,2 

(^VQ(0)^'-^Q'W^(0)^(t)), (31) 

where indices /i, ji (u, u) numerate the sublattices; = r il ~r v (r, is a vector indicating 
the CuOi centers within unit cell): 

is the Fourier component of the spin-charge susceptibility. 

Above expressions (24) - (31) enable to uncover a series of the unconventional properties 
of the induced spin fluctuations. Firstly, we note that the dynamics of the induced spin fluc- 
tuations can not only reflect but be determined by the dynamics of the charge fluctuations. 
This effect displays itself distinctly in the spin subsystem transparency region far from the 
main spin excitations, when we can neglect the temporal dependence of the spin correlation 
function in the expressions (29)-(31). Then, the temporal dependence of the induced spin 
fluctuations will be entirely determined by that of the charge fluctuations. 

Let's draw attention to the unconventional contribution in the induced spin fluctuations 
connected with the special for the singlet-triplet magnetic system one-center averages like 

(S m ) ((V m )) or in the ^-representation 

(S$S-?) = 5<ff{S ) 
In accordance with the common expression (29) at Q — Q', 

(6§M8$-&)) = \ p %\ 2 (S) (SN^0)6N_^(t)), (33) 

Q 



(32) 
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so the appropriate term in the induced spin fluctuations is independent on q and its t- 
dependence will be entirely determined by that of the charge fluctuations. 

The induced spin fluctuations result in unconventional contribution A% to the static 
magnetic susceptibility determined by the correlation function 

(^o(O)^o(O)) ~jjJ2\ pS (Q)\ 2 >< Cs (O)§_q(O))(5N (O)SN_ (O)) : (34) 

Q 

including purely spin correlation functions determining both the uniform (Q = 0) and stag- 
gered (Q 7^ 0) purely spin susceptibility. In this connection, note an appearance of the 
anomalously large "induced" spin susceptibility as a result of the similar " antiferromag- 
netic" fluctuations both in purely spin and charge subsystems. Actually, this corresponds to 
an appearance of the hybrid ferrimagnetic fluctuations (see Fig. 10). Peculiar feature of 
the induced ferrimagnetic contribution A^ to the paramagnetic susceptibility is its anoma- 
lous behavior near the temperature T C o of the charge ordering in the local boson system, 
schematically pictured in Fig.9. Note the "pseudo-gap" behavior at T < T C o, revealing 
in a relatively sharp decrease of the Ax with decreasing temperature below Too- So, the 
temperature Tco of the charge ordering will determine an universal temperature behavior 
of the induced ferrimagnetic contribution A% to the paramagnetic susceptibility typical for 
the critical regions. 

An interesting peculiarity of the induced spin fluctuations in the "bi-plane" cuprates like 
Y BaiCuj,0§+ x is an appearance of the ^-dependence of the correlation function (31) even in 
the absence of the purely spin exchange coupling between the nearest Cw02-planes. Really, 
the purely spin exchange interaction between the nearest C-uC^-planes is relatively weak 
that implies the possibility to neglect the terms with \i ^ ji' in (31) and so, in practice, the 

spin correlation functions like (S^^S^' ) for the PJT centers system do not depend on q z . 
On the other hand, the long range inter-plane Coulomb coupling and the appropriate boson 
transfer are relatively large that results in an essential Q z -dependence of the charge correla- 

tion functions like (SN^dN _^,) . Namely, this effect for the bi-plane cuprates is accompanied 
by a separation of the charge modes to the acoustical and optical ones in optimally doping 
compounds. So, these modes will be revealed in the induced spin fluctuations determining 

the corresponding (^-dependence of the correlation function like {SS^(0)SS-^(t)) ~ sin 2 (%^) 
or cos 2 (%p) (d is the distance between the C11O2 planes). 

For the description of the q, a> dependencies of the induced spin fluctuations it may be 
used approximations like the known Ornstein-Zernike model applicable to the systems with 
well developed fluctuations. Within a semiempirical approach we represent correlation func- 
tion (29) in the integral form as 

(<^0)<^_^)> = / ^(Q)(^-Q(0)l^(t)^Q(0)5iV_ (5 (t))dg, (35) 

where for the spin and charge correlation functions we can make use of the Lorentz or Gauss 
model form 

, 1 
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Here £ sp , ^ c h (a sp 2 , a c/l 2 ) are the correlation lengths for the spin and charge fluctuations, 
respectively. The resultant g-dependence of the induced spin fluctuations can have a remark- 
able variety of forms exhibiting the features of both spin and charge fluctuations. Fig. 11 
shows some typical curves modelling in the Gauss approximation the well known "four- 
peak" structure of the dynamical spin susceptibilities revealed by the inelastic magnetic 
neutron scattering for the cuprates like La2- x Sr x CuO i |TB[ . Incidentally, we conjectured 
an occurence of the antiferromagnetic singlet-triplet fluctuations with the specific wave vec- 
tor q sp = (1, 1) and an appearance of the the substitution driven charge or charge- vibronic 
incommensurate modes with q c h = (±2e, 0) or (0, ±2e) compatible with the stripe like su- 
perstructures with e 



x 
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5 The spin statics and dynamics in cuprates. 

An analysis of the various experimental data on magnetic susceptibility, neutron scattering, 
magnetic resonance (ESR, /xSR, NQR, NMR), Raman effect in the different cuprate HTSC's 
shows an absence of the long range magnetic order along with an existence of the well devel- 
oped spin fluctuations and sufficiently strong exchange interactions. As a rule, both static 
and dynamical spin susceptibility probed by different methods exhibit more or less distinctly 
a critical or "gap" behavior near the superconducting transition temperature T c and/or at a 
noticeably more higher temperature T* usually called as the "pseudo-gap" opening temper- 
ature. This phenomenon evidences the strong spin-charge intermode coupling. However, in 
some cases the superconducting transition doesn't result in either observable changes of the 
spin statics and/or dynamics. 

A complicated unconventional behavior of the spin system in cuprates can be in natural 
way explained at least qualitatively within the developed singlet-triplet PJT centers model. 
Below, within a framework of this model we'll shortly discuss a number of the most significant 
and hotly debated experimental results. 

Static spin susceptibility. As an example we'll consider the temperature and con- 



centration behavior of the static spin susceptibility xs m L(i2- x Sr x Cu0 4: [16], [LT]] which is 



a simplest model system among the HTSC's. The Fig. 12 schematically shows the typical 



Xs(T) dependences for a number of the 214-systems based on the experimental data [17] 
for x in the range 0.06 -v- 0.26. These can be reasonably explained in the framework of the 
model approach to the cuprates as to the singlet-triplet systems with strong fluctuations 
of the singlet-triplet separation A$t at the positive sign and relatively small magnitude of 
the corresponding average value A$t > or, in other words, at the predominantly singlet 
bare ground state of the PJT centers. Character of the temperature dependences Xs{T) 
for different x distinctly evidences the rise of these singlet-triplet fluctuations with the in- 
crease in concentration x accompanied by the possible decrease of the appropriate Ast 
value. Moreover, this explains a decrease in the T c value at x > x opt « 0.15 due to more 
and more effective spin reduction of the local boson transfer integral. In terms of the above 
developed singlet-triplet model the unusual temperature and concentration dependences of 
the paramagnetic susceptibility for L(i2- x Sr x CuO± in a wide range of temperatures T > T c 
and concentrations 0.06 4- 0.26 can be easily explained by the occurence of the SS — VV 
crossover regime for the singlet-triplet fluctuations. Observed decrease of xs with lowering 
the temperature below T max (x) can be coupled with the opening of the spin gap (see the 



inset to Fig. 5) in accordance with the conjecture of the authors |16|] that "the exceeding 
decrease of Xf(T) at low temperatures will be due to a development of some kind of singlet 
state". Note that an appearance of the appreciably large Curie term in the temperature 



14 



or 



behavior of the spin susceptibility xs (for example in La2- x Sr x CuO^ above x > 0.2 ]TJ 
in Tl 2 Ba2Cu0 6+x |Tj|) can be linked not with paramagnetic impurities but with an appre- 
ciable weight of fluctuations resulting in the paramagnetic "S'-type" contribution of the PJT 
centers. 

The minimal spin reduction of the local boson transfer integral essential for the opti- 



mized systems |T4| presumes a stability of the spin state of the PJT centers with regard to 
charge fluctuations. In this connection, one should pay attention to the temperature stable 
systems with the temperature independent normal state spin susceptibility as to the possible 
optimized systems with maximal T c . In such systems change in the temperature or a heat 
fluctuation doesn't accompanied by an appreciable change of the local spin state. Some 



authors |TI| consider condition ^ « as essential signature of the optimized systems like 
Y 'Ba2Cu3Oe.de- Note the interesting analogy with the vibronic reduction of the local boson 
transfer integral and condition of the minimal isotope-effect and baric coefficient ^ for the 



optimized systems ]T4| 



The Knight shift. The spin part of the Knight shift 11 K in the framework of the singlet- 



triplet model unlike the usually applied single-component spin model p0| includes two con- 
tributions connected with the usual spin susceptibility xs and unconventional susceptibility 
Xv- So, for the copper nuclei 63,65 Cti the taking account only of the direct hyperfine coupling 
gives 

63 ' 65 K sp = 7 r^—:A Cu ( X s + Xv), 

where Aq u is the direct hyperfine coupling parameter, fig is the Bohr magneton, and N is 
the Avogadro's number. On the other hand, for the 17 nuclei the taking account of the 
direct hyperfine coupling with the e u -orbital gives 

where Aq is the effective hyperfine coupling parameter. So, a comparison of the temperature 
dependences of the 63 > 65 K and 17 K allows to examine the singlet-triplet model and separate 
the S- and V^-contributions. As an illustration note the observed discrepancy between the 
temperature dependences of the 63 > 65 K and 17 K in La2- x Sr x CuO^ |H] that can be considered 



as a striking evidence in favour of the two-spin-component singlet-triplet model with the S 
and V contributions to the Knight shift. Note, that to a certain extent taking account of 
S and V terms is equivalent to that of localized copper spin and delocalized oxygen hole 
spin contributions sometimes enlisted for detailed analysis of hyperfine coupling and nuclear 
resonance in cuprates [p^j . 

Spin relaxation. Experimental studying the low-frequency spin dynamics in cuprates 
is mainly provided by the measurements of the NQR, NMR spin relaxation rates for different 
nuclei and the ESR spin relaxation rates for the rare earth ions (first of all for Gd 3+ ). 

Spin relaxation rates are mainly determined by the long range and long lived fluctuations 
that implies their critical (or near-critical) behavior near both the uniform phase transition 
temperatures like T/v, Tco, Tbs, T c and percolative phase transitions. Especial attention 
should be paid to the phase separation regime when a phase transition within single phase 
is obligatory accompanied by the long range and long lived fluctuations of the coexisting 
phases. 

A specific for the local boson model phase transition at T = Tco with the charge ordering 
at Q = (1, 1) should be revealed in a critical behavior of the spin-lattice relaxation rate Tf 
only for the copper nuclei 63,65 Cu, but not for the nuclei like l7 0, which crystallographic 
position results in a cancelation of the " antiferromagnetic" like fluctuations contribution. 
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Indeed, an expected critical behavior of T x 1 (or {T{T) l ) for the copper nuclei 63 ' 65 Cii 
has been observed in La 2 _ x Sr x CuO± (Tco ~ 50 K for x = 0.13 and x = 0.18 P5|), in 
YBa 2 Cu 3 6m (T co « 150 AT 0), in optimized YBa 2 Cu 3 6 . 96 , where T co ~ T BS « T c 



20f1 . Appropriate temperature dependences for the ^(TiT)^ 1 within the above developed 
induced spin fluctuations model in accordance with expressions like (29) are determined by 
the temperature dependences of the charge fluctuations (SNaSN a)u>-*q an d look like those 
shown in Fig. 9. 

A relatively low value of the spin-lattice relaxation rates for the 17 0, 89 Y, 139 La and their 
temperature dependences can be qualitatively explained within a singlet-triplet model with 
assumptions used in the case of static spin susceptibility. 

Inelastic neutron scattering. Magnetic inelastic neutron scattering is a powerful tool 
for studying the detailed u, ^-dependences of spin fluctuations. An appearance of new phase 
of the PJT centers is accompanied by a dramatic modification in magnetic excitations with 
large broadening in wave vector and substantial redistribution of spectral weight in frequency 
p4j| . Magnetic fluctuations extend to a wide spectral range up to energies ~ 0.3 eV with 



spectral weight removed from the appropriate magnetic Bragg peak observed for the parent 
antiferromagnetic compounds. To a certain extent the magnetic response for the basic 
model HTSC's La 2 - x Sr x Cu04 and YBo 2 Cu 3 Oq +x can be represented as a superposition 
of a strong broad background due to antiferromagnetic fluctuations with Q = (1, 1) and 
the relatively weak and sharp low frequency generally incommensurate features peaked near 
lOmeV in 214 system and near 30 4- AOmeV in 123 system. A structure of the basic 
" antiferromagnetic" background changes with concentration x exhibiting a phase separation 
regime with a suppression of the parent phase contribution and a corresponding rise of 
the PJT centers phase singlet-triplet contribution. Antiferromagnetic correlation length for 
the former decreases with doping from relatively large values for slightly doped compounds 
to minimal values £ ~ RcuCu, whereas that of the latter remains minimal indicating a 
permanently short-range character of the antiferromagnetic singlet-triplet fluctuations. As 
seen from the experimental data [24] a normalized intensity of the magnetic inelastic neutron 
scattering for the PJT centers phase is in common sufficiently weaker than that for the 
parent antiferromagnetic phase at least at T < 300 K . This is entirely compatible with the 
qualitative conclusions of the singlet-triplet model with predominant lty singlet ground state 
of the PJT centers. Unfortunately, the experimental data for the YBo 2 Cu 3 Oq +x are limited 
to relatively low frequencies (hu ~ 120 meV) that doesn't permit to correctly examine the 
singlet-triplet contribution. Nevertheless, the available magnetic INS-data for various doping 
level agree with a superposition model within phase separation regime and can be effectively 
treated in framework of the two Lorentz oscillators model [|T5]]. 

Phase separation and magnetic properties of cuprates. Phase separation both 
chemical and physical should be taken into account for a correct interpretation of the avail- 
able experimental data. Properly speaking, within a phase separation regime we deal with 
contributions of two (or even more) phases with varying relative volume. In many cases, 
there are some difficulties with a distinct separation between single-phase and multi-phase 
effects. In particular, it concerns the phase transition phenomena accompanied by the vary- 
ing the relative phase volume. Moreover, a phase separation regime can exhibit a specific 
phase transition of the percolative origin. 

Below, we list shortly some experimental results of the static and dynamic magnetic 
measurements which more or less distinctly indicate an occurence of the phase separation 
regime in either cuprates. 

As one of the recent remarkable experimental indications to the nucleation of the PJT 
centers phase in the cuprates note the zero field copper NMR data in Y\_ X C 'a x Ba 2 Cu 3 0% 



16 



by P. Mendels et al. ||25|| . The nonisovalent substitution in the antiferromagnetic state was 
accompanied by the anomalous decrease in the concentration of the NMR resonating copper 
nuclei: every Ca 2+ -ion excluded from the NMR about 50 copper ions (!), that could be 
connected with their disproportionation within the new phase "droplets". A decrease in 
the volume of the antiferromagnetic phase in YBaiCujQ^+x with increasing x results in a 
suppression of the Bragg peaks intensity in the inelastic neutron scattering spectra [26]. 

An occurence of two characteristic temperatures, namelly, that of the charge ordering 
{Too) and appearance of the local superconducting order (Tbs) gives rise for the two- 
dimensional regions of the PJT centers phase within the CtiCVplanes. It can be revealed in 
many effects even far from the percolative transition to the superconducting state x < 0.06 
for Ld2-xSr x Cu04 and x < 0.45 for Y Ba2Cu^O^ +x or within a formally antiferromagnetic di- 
electric state. So, a well known anomalous temperature behavior of the Bragg peak intensity 
and spin-echo relaxation rate in Y Ba2Cu%0§ +x PB], |27J ("reentrant transition"), schemati- 
cally shown in Fig. 13, can be explained by a decrease in a volume of the antiferromagnetic 
regions with decreasing the temperature below T ~ 30 K (Tag?) ■ A similar effect is revealed 
by an anomalous temperature dependence of the Zeeman splitting in the 139 La NQR spectra 
for La2~ x Sr x CuO^ at T < 30 K ]p8[] . A percolative nature of the superconducting transition 
in La2- x Sr x Cu04 at x ~ 0.06 reveals in anomalous rise of the low temperature 139 La spin 
relaxation rates [29[ . 



Interestingly, that strong isotop substitution effect in the spin susceptibility and the Neel 
temperature observed for La,2- x Sr x CuO^ [0 can be linked with an influence of the isotope 
substitution 16 — > 18 O on the relative phase volume. 



6 Conclusions. 

A suggested singlet-triplet model considers the spin (magnetic) properties of the cuprates as 
an integral part of their multi-mode behavior characteristic for the systems with suficiently 
strong inter-mode coupling. Incidentally, it is rather simplified approach. In particular, 
we did not touch seriously upon the problem of the vibronic states and appropriate local 
structural Jahn- Teller modes. Actually, in cuprates we should deal with hybrid spin-charge- 
vibronic fluctuations. 

Another important problem typical for the high-T c cuprates is connected with their quasi- 
two- dimensionality and its influence on the spin arrangement. 

A supplementary important issue is coupled with the charge inhomogeneity and phase 
separation. These phenomena should be taken into account for an interpretation of the 
temperature and concentration dependencies of all both static and dynamical spin properties, 
especially, within the range of the critical behavior of either subsystem accompanied by the 
change of the relative weights and/or volumes of the coexisting phases. 

Despite the used simplifications the singlet-triplet model enables to make some qualitative 
and semi-quantitative conclusions compatible with the observed experimental data. First of 
all, it concerns an occurence of the singlet-triplet antiferromagnetic fluctuations of various S 
and V types resulting in various temperature dependence of different quantities describing 
the spin static and dynamic properties. Let note, that unlike the widely popularized nearly 
antiferromagnetic Fermi-liquid model |2D[ the suggested approach naturally enables to take 
account both of d- and p-holes in conditions of their comparable role. 

As a whole, singlet-triplet model can provide a promising foundation for the qualitative 
and semi-quantitative analysis of the static and dynamical spin properties of cuprates includ- 
ing paramagnetic susceptibility, nuclear resonance and inelastic magnetic neutron scattering. 
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Moreover, the further elaboration of the singlet-triplet model represents a perspective task 
not only for the HTSC-problem. 
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Figure captions. 



Fig.l. The structure of the lower energy levels for the PJT center. The four types 
of distorsions of the CuO^ clusters corresponding to the four minima of the ground state 
adiabatic potential are schematically depicted in the insets for the B\ g and B 2g symmetry 
types, respectively. 

Fig. 2. The electron density distribution for the 1,3 E U term of the b\ g e u (a) configuration 
(left) and the bi g e u (-K) one (right) compatible with the B\ g type distorsion of the CuO^ 
clusters. The dark filling corresponds to the e u orbitals and light filling does to the b lg ones. 
The appropriate vibronic states are ascribed to the a z = ±1 pseudo-spin states. 

Fig.3. Illustration to possible angular (a) and collinear (b,c) spin S, V configurations for 
the singlet-triplet magnet. 

Fig.4. The orientations of the copper and oxygen spins for the different collinear phases 
(qualitatively). The distinguishing tone filling corresponds to the possible different orbital 
configurations. 

Fig.5. The temperature dependence of the transversal (solid curve) and longitudinal 
(dotted curve) spin susceptibilities xs,v an d order parameters for the SS phase (A$t = 
O.leV, zlW = OAeV, zl<® = -O.SeV, zl& = 0.2eV). The inset shows an appropriate 
simplified spin level scheme (T = 0). 

Fig.6. The temperature dependence of the transversal (solid curve) and longitudinal 
(dotted curve) spin susceptibilities xs,v an d order parameters for the VV phase (A S t = 
O.leV, zlW = 0.2eV, zl& = -0.3eV, zl& = OAeV). The inset shows an appropriate 
simplified spin level scheme (T = 0) with the spin gap. 

Fig. 7. (a) The temperature dependence of order parameters for noncollinear phase (A = 
O.leV, zl {1) = 0.2eV, zl {2) = 0.3eV, zl^ = -O.leV); (b) The temperature dependence of 
"intra-site" Q A a{T) = ©rb(T) and "inter-site" ®ab{T) angles (see text). 

Fig. 8. The temperature dependence of longitudinal (solid curve) and transversal (dotted 
curve) magnet susceptibilities Xs,v f° r noncollinear phase (A = O.leV, zl^ = 0.2eV, zl^ = 
0.3eV, zlW = -O.leV). 

Fig.9. (a): A model (T, z/)-phase diagram for the quantum lattice bose-gas (N B ~ 0.5) 
with v being a parameter (in particular N B or doping content!) determining a certain change 
of the boson-boson repulsion V and the boson transfer integral t (see the inset); (b, c, d) show 
a qualitative T-dependence of the charge fluctuations and the superconducting order param- 
eter fluctuations for three regimes marked by the arrows in (a): V > t (weak screening), 
V m t (optimal screening) and V < t (over-screening). The standard abbreviations denote 
the non-ordered or metallic phase (NO), the charge order phase (CO), the bose-superfluid 
(superconducting) phase (BS), and mixed superconducting) phase (BS+CO). 

Fig. 10. An illustration to the induced ferrimagnetic ordering within the singlet-triplet 
model. Different filling for the CuO^ centers indicates different charge (boson) density. For 
the sake of the close approach to the real cuprates like La 2 - x Sr x Cu04 we indicate also the 
buckling of the Cu0 2 planes (by different filling of the appropriate oxygen ions) and possible 
different sizes of the CuO^ clusters with different charge density. 

Fig. 11. The model "four-peak" structure of the dynamical spin susceptibility for the 
cuprates like La 2 - x Sr x CuO^ determined by induced spin fluctuations and revealed by the 
inelastic magnetic neutron scattering. Incidentally, it was conjectured an occurrence of the 
antiferromagnetic singlet-triplet fluctuations with the specific wave vector q sp — (1, 1) and an 
appearance of the the substitution driven charge or charge-vibronic incommensurate modes 
with q ch = (±2e, 0) or (0, ±2e) compatible with the stripe like superstructures with e ~ x. 
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Fig. 12. The typical Xs(T) dependences for a number of the 214-systems based on the 
experimental data |TB|, |T7J. The curves 1-6 correspond to x increasing in the range 0. 06-^0.26. 
Arrows indicate T max . The inset shows a qualitative energy spectrum for the PJT centers 
with taking account of the well developed fluctuations. 

Fig. 13. Anomalous temperature behavior of the Bragg peak intensity and spin-echo 
relaxation rate in Y Ba2- x Cu^OQ_^) [27j near T = 30 K (schematically). 



21 



This figure "figl.gif" is available in "gif" format from: 



http://arXiv.org/ps/cond-mat/971 1 1 16vl 



This figure "fig2.gif" is available in "gif" format from: 



http://arXiv.org/ps/cond-mat/971 1 1 16vl 



This figure "fig3.gif" is available in "gif" format from: 



http://arXiv.org/ps/cond-mat/971 1 1 16vl 



This figure "fig4.gif" is available in "gif" format from: 



http://arXiv.org/ps/cond-mat/971 1 1 16vl 



This figure "fig5.gif" is available in "gif" format from: 



http://arXiv.org/ps/cond-mat/971 1 1 16vl 



This figure "fig6.gif" is available in "gif" format from: 



http://arXiv.org/ps/cond-mat/971 1 1 16vl 



This figure "fig7.gif" is available in "gif" format from: 



http://arXiv.org/ps/cond-mat/971 1 1 16vl 



This figure "fig8.gif" is available in "gif" format from: 



http://arXiv.org/ps/cond-mat/971 1 1 16vl 



This figure "fig9.gif" is available in "gif" format from: 



http://arXiv.org/ps/cond-mat/971 1 1 16vl 



This figure "figl0.gif" is available in "gif" format from: 



http://arXiv.org/ps/cond-mat/971 1 1 16vl 



This figure "figll.gif" is available in "gif" format from: 



http://arXiv.org/ps/cond-mat/971 1 1 16vl 



This figure "figl2.gif" is available in "gif" format from: 



http://arXiv.org/ps/cond-mat/971 1 1 16vl 



This figure "figl3.gif" is available in "gif" format from: 



http://arXiv.org/ps/cond-mat/971 1 1 16vl 



